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QCD lattice simulations yield hadron masses as functions of the quark masses. From the gradients 
of the hadron masses the sigma terms can then be determined. We consider here dynamical 2+1 
flavour simulations, in which we start from a point of the flavour symmetric line and then keep 
the singlet or average quark mass fixed as we approach the physical point. This leads to highly 
constrained fits for hadron masses in a multiplet. The gradient of this path for a hadron mass then 
gives a relation between the light and strange sigma terms. A further relation can be found from 
the change in the singlet quark mass along the flavour symmetric line. This enables light and 
strange sigma terms to be estimated for the baryon octet. 
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1. Introduction 

In this talk we shall describe a method for the determination of the hyperon sigma terms based 
on the results of [|I|] to which we refer to for more details including numerical results. 

(H) (H) 

Sigma terms, , a} are defined as that part of the mass of the hadron (for example the 
nucleon) coming from the vacuum connected expectation value of the up (u) down (d) and strange 
(s) quark mass terms in the QCD Hamiltonian, 

a/^) =m'}{H\ {uu + ddf \H) , ai^^ = < I {ss f \H) , (1.1) 

where we have taken the u and d quarks to be mass degenerate, niu = nid = mi. (The superscript " 



denotes a renormalised quantity.) Eq. (1.1) is usually written (in particular for the nucleon) as 



{N) _ m1{N\{uu + dd-2-ssY\N) ^j^y _ 2{N\{^sf\N) 

^' ~ l-yW'< ' ^ ~ {N\{uu + ddY\N) 



(i.e. we consider y^^^" rather than ai^'). The simplest calculation, (which we will discuss in more 
detail later) uses first order in SU (3) flavour symmetry (octet) breaking to give 

(N) mf M-w+Mi,-2Mn 26 
a}^> = ' " ,,,, — ^MeV, (1.3) 



and 



where / ni^ is the ratio of the strange to light quark masses, which using the leading order PC AC 
formula for this ratio gives m^/mf = {2MI-MI)/MI ~ 25. The Zweig rule, {N\{ssy\N) ~ 
would then give a/^^ ~ 26MeV, ai^' ~ OMeV while any non-zero strangeness content, y^^^" > 
would increase this value of Cj/^\ ai^^ (and indeed due to the large coefficient, ai^^ quite rapidly). 

Computing the sigma terms from lattice QCD has a long history from quenched to 2 flavour 
and more recently 2+1 flavour simulations. In general more recent results tend to give lower 
values than earlier determinations. 

In this talk, we shall investigate this simple picture as described above and in particular test 
the linearity assumption of SU (3) flavour symmetry breaking. 



2. Flavour symmetry expansions 

Lattice simulations start at some point in the (mj,m^) plane and then approach the physical 
point (m^*,m^*) along some path. (In future we shall denote the physical point with a *.) As we 
shall be considering flavour symmetry breaking then we shall start here at a point on the flavour 
symmetric line ml = mf and then consider the path keeping the average quark mass constant, 
m = const.. The SU{3) flavour group (and quark permutation symmetry) then restricts the quark 
mass polynomials that are allowed, giving for the baryon octet 

Mh =Mo{m)+CH5mi + 0{5mi), (2.1) 
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with cn = 3Ai, ca = 3A2, = —SAj, 



-3(Ai -A2) and 



5mi=mi — m, m=4(2m;+m,s 



(2.2) 



So to linear order in the quark mass, we only have two unknowns, Ai, A2 (rather than four). A sim- 
ilar situation also holds for the pseudoscalar and vector octets (one unknown) and baryon decuplet 
(also one unknown). This highly constrains the numerical fits. 

Permutation invariant functions of the masses, Xs, (or 'centre of mass' of the multiplet) have 
no linear dependence on the quark mass. For example for the baryon octet we have 



Xn = ^{Mn+Mz + M2)= Mo (m) + 0{5mj 



(2.3) 



(The corresponding result for the pseudoscalar octet is given later in eq. ( 3.10| ).) 

Furthermore expanding about a specific fixed point, m/ = = mo on the flavour symmetric 
line and allowing m to vary, we then have 



Mo{m) = Mo (mo) +Mo(mo)(m - mo) + 0{{m- mo) 



(2.4) 



We will see that Ai, A2 determine all the non-singlet sigma terms and M'{mo) the singlet sigma 
terms. 

As an example of the quark mass expansion from a point on the flavour symmetric line in 
Fig. [I| we plot the baryon octet Mh/Xn for H = N,A,'L,E against M^/X^ together with a linear fit. 
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Figure 1: Mh /Xn (H = N, A, E, S) against M'^/X^ for initial point on the flavour symmetric line given by 
Kb = 0.12090, left panel and fCo = 0.12092, right panel. The 32^ x 64 lattices ai-e filled circles, while the 
24^ X 48 lattices are open triangles. Also shown is the combined fit of eq. (3.11) (the dashed lines) to the 



32 X 64 lattice data. The fit results are the open circles, while the experimental points are the (red) stars. 



eq. (lA) and implicitly eq. ( 3.1C ) using 2+1 0{a) improved clover fermions at j3 = 5.50, using two 
starting values for the quark mass on the flavour symmetric line. All the points have been arranged 
in the simulation to have constant m. We see that a linear fit provides a good description of the 
numerical data from the symmetric point (where Mji ~ X * = 410.9MeV) down to the physical 
pion mass. 
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3. (Hyperon) o equations 
3.1 Renormalisation 

For Wilson (clover) fermions under renormalisation the singlet and non-singlet pieces of the 
quark mass renormalise differently. We have 



m\ = Z'' [niq + az\{2mi + m,)] , ttz 



(3.1) 



In the action the term Y,y'^qQQ = Y^q^qi^llY i-^- ^ renormalisation group invariant or RGI quan- 
tity. Upon writing this in a matrix form and inverting gives 



qq 



az 1 



{uu + dd + ss) 



(3.2) 



so for az / then there is always mixing between bare operators. Useful quark combinations are 
the octet^ and singlet combinations, namely 



{uu + ddf -2{^sf 



^ {{uu + dd)-2{ss)) 
{m + ddf + [ssf = 2Ns(^^^az) ^ ^"^^ ^ ^''^^ 



(3.3) 
(3.4) 



3.2 a equations 



Scalar matrix elements can be determined from the gradient of the hadron mass (with respect 
to the quark mass) by using the Feynman-Hellman theorem which is true for both bare and renor- 
malised quantities. So if we take the derivative with respect to the bare quark mass we get the bare 
qq matrix element, 



H 



dmi 



{H\{uu + dd)\H) 



dM, 



H 



{H\ss\H) . 



(3.5) 



Multiplying the renormalised quark mass, eq. (3.1), together with eqs. (3.4) (or more generally 
with eq. (^^)) we can find RGI combinations (i.e. a form where the renormalisation constant 
cancels). In particular we find 



2rol 



3r 



l+2r 
3r 

T+2r 



{l + az)moCH 
moM,)(mo), 



(3.6) 
(3.7) 



where r is the ratio of quark masses r = mf/m"^. The two simultaneous equations, which can be 
easily solved, give 



l+2r 
1 

l+2r 



[(1 + az)moCH + 2moMQ{mo)] 
[-(1 + az)moCH + moMo(mo)] 



(3.8) 



'Eq. (y)also leads to eq. (|l.3|) as from section g^we have {N\{uu + dd -IssflN) = c^/Z'^'^ = 3Ai /Z^'^ . Together 
withMs +Mx — 2Mf^ = — 9Ai5m/ = 3Ai{in^ — mf)/Z'^^ this gives eq. (L3l. An alternative mass combination that also 
picks out the Ai coefficient is Ms — M/^ = — 3Ai 5m/. 
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(H) (H) 

We see that the smallness of in comparison to a) is certainly guaranteed by the presence of 
an additional r in its numerator. These coefficients are also sufficient to determine yW"^ as can be 



seen from eq. ([3.8|), 



y(H)R ^ 2 -(1 + az)mQCH + moM'Q{mo) 
(1 + az)moC// + 2moM,',(nio) ' 

It is seen that y^^'^" only depends on gradients and not on the physical point. 

It is convenient to normalise the coefficients by Xj^ so we now need to find the coefficients 
(1 + az)mQCH /Xa?(wio) and moMQ^mo) /XA?(mo) and also to extrapolate to the point where the quark 
mass ratio takes its physical value, i.e. r = r*. 

3.3 Determination of the coefficients 

The coefficients can be determined by considering gradients with respect to a physical quantity. 
As in eq. ( ^?T| ) we also have a similar expansion for the pseudoscalar octet, 

Ml = mIj^ + 2a5mi + 0{8mj) , (3. 10) 

(together with M|- = Mq^^ - admj + 0{dmf), M^^ = M^^j^. -Aa5mi + 0{dmj)). Analogously to 



eq. ( |2.3| ) we can define a flavour singlet quantity X| = \{2M\ +M|) = Mq^j. + 0(5mj) However, 
as well as eq. (^!1|), we have the additional constraint from PCAC = IBf^ml (together with 
Ml =Bf,(m^ + mf)). If we now consider an expansion in the pion mass then eliminating 5mi 



between eq. (Ejl) and eq. (3.10) gives 



= (l - [(1 + az)mo"^]) + [(1 + az)mo'^] ^ , (3.1 1) 

An \ Am / An Aj^ 

from the point on the symmetric line mo = m. Thus if we plot Mh/Xn versus M^/X^ (holding the 
singlet quark mass, m constant) then the gradient immediately yields (1 + az)moCH /X^. The only 
assumption is that the 'fan' plot splittings remain linear in dmi down to the physical point. In Fig. |] 
we show this plot leading to a results for (1 + az)moCH /X^ for K"o = 0. 12090, 0. 12092. 



Furthermore on the flavour symmetric line eliminating (m — mo) between eqs. (2^) and the 
corresponding one for M^(m) gives 

X/v(m) _ ^ moM^(mo) A ^ ^ moM^(mo) j X|(m) 



Again in a plot of XN{ln)/XN{mo) versus X|(m)/X|(mo) the gradient immediately gives the re- 
quired ratio moMQ(mo)/XA?(mo). In Fig. ^ we plot XN{rn)/XN{mo) versus X|(m)/X|(mo). From 



eq. (3.12) the gradient gives the required number. 



Finally the quark mass ratio, r, must be estimated. In the right panel of Fig. ^ we plot (2M|- - 



M„)/X^ versus MyX^. From eq. (3.1C) we have 



2Ml-Ml ML Ml 

— ^ — ^ = 3_^_2^. (3.13) 

An Aj^ Af^ 

As in section ^, we see that for constant m the data points lie on a straight line (i.e. there is an 
absence of significant non-linearity). Together with PCAC, this gives the A;-axis is proportional to 
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Figure 2: Left ploV.XN{m) /XN{mo) versus X^{m)/X^{mo) along the flavour symmetric line, together with 
flie linear fit from eq. ( [3.12| ); Right plot: {IM^-MD/X^ versus mI/X^ for Kq ^ 0.12090 (left panel) and 



Kb — 0.12092 (right panel). The 32^ x 64 volume results are given by the filled symbols, while the 24^ x 48 



volume results are shown using empty triangles. The fit is given in eq. (3.13). Experimental points are 
denoted by (red) stars. 



m'l while the y-axis is proportional to and thus the ratio gives r. Taking our physical scale to be 
defined from M^/X^j* (i.e. from the x-axes of Fig. ^) gives 1 /r*. 

What can we say about corrections to the linear terms? The simple linear fit describes the 
data well, from the symmetric point to our lightest pion mass, both along the rn = const, line and 
the flavour symmetric line. To see the possible influence of curvature we compare linear fits with 
quadratic fits as discussed in the Appendix of [|T]]. These will be used in section^ for the estimate 
of possible systematic effects. 



4. Results and Conclusions 



We can now determine y'^^^" and oj^K Gs^\ The scale is taken asX^ = l.lSOlGeV. We shall 
only discuss here the general details of the results; the numerical values are given in 



From eq. ( |3^ we can find an indication of the magnitude of a/^' as approximately 



a/"''*~22.4 + ^j^MeV>22.4MeV. (4.1) 

The last inequahty follows as obviously Os > 0. Indeed this shows that a non-zero Os > 
can only add a few MeV to this result. 

These results are illustrated in the left plot of Fig. ^ for y^^^'^* and in the right plots of^^"' , 
Os^^* both against H = N, A, L and E. While the data for Kq = 0.12090 is more complete than 
for Kfl = 0. 12092 (cf. the plots in Fig. |T]) and demonstrates linear behaviour, as the path starting at 
Kq = 0. 12092 is closer to the physical point (cf. Fig. ^ we shall use these values as our final values. 

In conclusion we have found that keeping the average quark mass constant gives very linear 
'fan' plots from the flavour symmetric point down to the physical point. This implies that an 
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Figure 3: yy"'"* (left plot) and a/ ' and a, ' (right plots) both at the physical point for H = N,A,'L,Z 



expansion in the quark mass from the flavour symmetric point will give information about the 
physical point. In this talk we have apphed this to estimating the sigma terms (both light and 
strange) of the nucleon octet. There has been no use of a chiral perturbation expansion (indeed this 
is an opposite expansion to the one used here, expanding about zero quark mass). 

Note that expansions about the SU (3) flavour line require consistency between many QCD 
observables, here for example not only for the baryon octet under consideration here, but also for 
the pseudoscalar octet, and PCAC and the ratio of the light to strange quark mass. 

Of course there are several more avenues to investigate. Our approach here has been to em- 
phasise linearity at the expense (presently) of reaching exactly the physical point. This can be 
addressed by interpolating between a small set of constant m lines about the physical point. Ad- 
ditionally the use of partial quenching will also help to get closer to the physical pion mass. With 
more data, a systematic investigation of quadratic quark mass terms in the flavour expansion should 
be considered, to reduce the systematic errors. 
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